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Abstract 


The pass£if;e of debris frotn a high altitude (>400 km) nuclear burst over the 
ionospheric pk.sma is found to be capable of exciting large amplitude whistler 
waves which caft act to structure a collisiohless shock. This instability will occur 
in the loss cone exits of the nuclear debris bubble^ and the accelerated ambient 
ions will free-stream along the magnetic field lines into the magnetosphere. Us- 
ing Starfish-like parameters and accounting for plasma diffusion and thermalization 
of the propagating plasma mass, it is fotmd that synchronous orbit plasma fluxes of 
high temperature electrons (near 10 keV) will be significantly greater than those 
encountered during magnetospheric substorms. These fluxes will last for suffi- 
ciently long periods of time so as to charge immersed bodies to high potentials and 
arc discharges to take place. Synchronous orbit satellites expecting to operate in 
a high latitude, hi'^h altitude nuclear burst environment should be designed against 
this effect as well as the radiation encountered. 
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1. IM UOIM ( IION 


Turbulent coupUng between the ejected debris plasma and background air 
plasma of a high-altitude nuclear burst appears capable of bringing about electron 
streaming to high altitudes. Ih the presence of the compressed .nagnetic field 
Which can penetrate the debris bubble, plasma turbulence takes the form of large 
amplitude whistler waves. These waves can couple the background air plasma 
with the expanding debris plasma. Fot electrons are acquired and some have 
sufficiently large velocities to escape along the distended field lines. The nonlin- 
ear wave-resonant particle interactions produce anomalous resistivity whose sea e 
length determines the extent of magnetic field penetration into the bubble which, in 
turn determines the rate of escape of (i and plasma electrons. 

’>3uperalfve'nlc debris plasma can escape directly through loss cone exits in the 
debris bubble since, ih general, the cylindrical axis of the bomb casing is no 
initially aligned with. the direction of the geomagnetic field. The more perpendic- 
ular the cylindrical axis is to the field lines, the greater the number of such 
escape particles. Their superalfvenic velocities suggest the formation of paralle 
collisionless shock waves ' geomagnetic field). Studies reveal 

that such shock fronts are structured by turbulent whistler modes which couple the 
incoming background air plasma to the -hocked debris plasma. Air plasma can 
therefore be picked up by the loss cone debri. and accelerated to high altitudes. 

The presence of a magnetic field has a significant effect on shock wave struc- 
ture. Gradients in the magnetic field give rise tc electron currents that 
ion acoustic Waves unstable and increase the effective collision frequency « 

(this dictates the penetration depth of the compressed magnetic field into the debris 
bubble, so that the rate of escape of debris and air electrons is profoundly affected). 
When propagation is perpendicular to the magnetic field, the magnetic field can 
inhibit the electrons from shorting out ion plasma oscillations for wavelengths long 
compared with the electron gyroradius‘^.;> and wavelengths short com- 
pared with the electron gyroradius. interactions between the ion beam mode 

and the electron Bernstein modes generate instabilities which are stabilized by 
electron heating, resonance broadening, or ion trapping. For oblique or parallel 
propagation, interactions of whistler waves with ion acoustic beam modes or 
Ion-cyclotron beam modes ^ are likely to be Important, and the existence of 

Whistlers depends upon the presence of a magnetic field. 

An instability found to be especially attractive as a colllsionlcss mechanism 
for pickup and heating of air electrons is the ion cyclotron beam modc-whlstler 
mode (current-free) instability that Golden^' found to be operative along the field 
lines and particle trajectories issuing directly from the loss cone exits. 
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The collisionless shock waves can be modeled as Mott-Smith layers of inter- 
penetrating unshocked (background air) and shocked (thermalized debris-air piston) 
flows, so that these layers are natural environhients for streaming instabilities. 
Recent linear dispersion investigations of parallel shock layers^ which form 
ahead of the debris plasma issuing from the loss cone exits reveal that, for a given 
Mach number > M* 2. 77 (M^ - = shock velocity, = Alfve'n 

speed in unshocked plasma), the shocked hot ion cyclotron beam mode can always 
drive unstable a particular whistler mode which, in the rest frame of the shock 
front, is stationary near the leading edge. An analysis of the shock interior 
reveals that the shock Mach number determines the portion of the shock thickness 
in which unstable whistlers are stationary hi the shock rest frame. For = 2.77, 
such modes may stand only at the leading edge, whereas for stronger shocks 
(M^ > 2. 77 ), they may stand at all points between the leading edge and some interior 
point which is dependent ort shock strength. For very strong shocks (M^ » 1), 
fully one -third of the shock thickness is filled with these modes, which can there- 
fore grow to large amplitude and couple the quiescent background air plasma to the 
expanding debris -air piston. 


2 . iNSniill.lTV I:\(;IT\T10N 

The number of such debris particles which enter a loss cohe exit depends 
critically on the mass and orientation of the bomb casing just prior to the burst. 

In the loss cone corridor, the debris plasma drives a shock wave. This shock is 
structured by unstable whistler waves which stand at its leading edge. To verify 
that these whistler modes can grow to sufficiently large amplitude to scatter incom- 
ing air ions (as viewed in the reference frame of the shock front) the following 
analysis is performed; During daylight burst conditions the density of oxygen ions 
is no + ~ 10 -10 , the larger value reflecting maximum sunspot activity. For an 
ambient field strength B ~ 0. 3 gauss, this corresponds to an Alfven speed 
Ca ~ 163 - 516 km/sec irt the ambient air plasma. 

For typical initial casing velocities of from 500 km/sec to 2000 km/sec, the 
initial Alfven Mach numbers are from 1. 96 to 12. 22. A summary of the spectrum 
of stationary unstable leading edge whistlers, based on the whistler dispersion 
relation, is given in l^lgure 1. It is seen there that a broad range of wave num- 
bers can grow to large amplitude during daytime (n ^ 10® cm‘^) bursts. To 
investigate whether these modes can achieve these large amplitudes in the loss 
cone exit, we must calculate the growth rate of these unstable waves. For simplic- 
ity we choose | 2 for which M^ 4. 12. indicating a realistic initial 

debrli^ vg/ city of u ()74 km/sec. In this case the calculated linear growth rate 
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Figure 1. Summary of Unstable Leading Edge Whistlers 
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for Icadtnii edge modes, from the whistler dispersion relation, is 7 
where rj is shock strength expressed as the ratio of shocked to unshocked densi- 
ties. For i “ 500 km and n = • 2 at the shock leading edge, we find 7 </u « lo, 
indicating 10 e-folding times pass as the shock traverses 500 km of the loss 
cone corridor. This 1 kl « unstable whistler mode can grow to suf- 

ficiently large amplitude after traveling 500 kn> along the loss cone corridor to 
pick up background air plasma. The air electrons are rapidly accelerated by 
ensuing ion-electron electrostatic instabilities. Since the calculation is based 
on the lower bound of it i® clear that 10 e-folding times irt 500 km is a con- 
servative estimate since Stronger shocks will have larger growth rates. 


i. 0i‘ 111 list I I I M S 'MTU MM.NKTOSiMllillH: M llSKIHM I l.l \i;> 

In the past few years there has been considerable concern over the phenomena 
of synchronous orbit satellites charging to high potentials as a result of magneto- 
spheric substorms. These substorms consist of the injection of high energy 
plasma from the earth's magnetotail into the region of synchronous orbit. Those 
portions of a satellite subject to the high energy plasma will charge to a potential 
several times the electron energy, while other portions of the satellite will remain 
at ground potential. Potentials near ground are maintained by photoelectron emis- 
sion from illuminated surfaces on the spacecraft or by contact with the ambient 
low energy plasma. 

During eclipse photoelectron emission disappears, and during a substorm the 
ambient low energy’ plasma flux is strongly dominated by the injected high energy 
plasma. The most damaging discharges probably occur between shadowed space- 
craft components influenced by substorm plasma and illuminated components at 
ground potential. M'hen the discharge passes through electrical circuitry between 
the components, damage can result. Electromagnetic interference can also result 
from surface discharges and considerable surface deterioration can be caused by 
arcing. 

The following discussion will be an assessment of the possible spacecraft 
charging effects which can result from the large scale transport of ionospheric 
plasma to synchronous orbit by a nuclear burst. The plasma instability just dis- 
cussed demonstrates a mechanism for structuring a collisionless shock wave. 

This mechanism will operate during a high altitude nuclear burst as the expanding 
nuclear debris passes over the stationary ionospheric plasma. Through the inter- 
action of large amplitude whistler mode waves, plasma will be picked up by the 
collisionless shock and accelerated into the magnetosphere. 
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Using nuclear burst parameters for a Starflsh-Uke (nominally 1. 5 megaton) 

burst, calculations^^ have conservatively estimated that near 10^^ 0^ ions 

would be carried by the shock when 1 percent of the total debris exits the loss 

cone. Depending on orientation, mass, and shape of the bomb casing, this number 

could be considerably higher. This number is reasonably estimated by the mass 

of debris Which exits the loss cone, for ion pickup slows the debris piston and 

eventually shuts off the pickup instability. 

The total propagation time for the plasma mass from the burst point just 

above the earth's surface to synchronous orbit is tens of seconds and more than 

sufficient to have thermalization of the complete plasma mass at near the ion 

17 

temperature as the mass slows and diffuses. 

27 4 

The estimated fiickup of 10 0 ioiis by the nuclear burst shock wave would 

have an energy of approximately 10 keV at reasonable shock velocities. The 

demands of plasma neutrality would quickly accelerate an equal number of electrons 

which would thef'malize with energies equal to or greater than the ion energy. 

At the loss Cone exits these 10^^ electrons with energies of 10 keV would be in a 

21 3 

cylindrical volume of approximately 2 y lO^^ cm , assuming 1 percent of the 

nuclear burst bubble as comprising the loss cone exit and the p;.asma pickup region 

being several hundred kilometers in extent. These numbers and energies translate 
° 15 2 

to an omnidirectional flux at the loss cones of * 3 v 10 electrons /cm -sec. 

The calculation of fluxes at higher altitudes than the burst altitude follows 

directly from the Liouville Theorem that a differential intensity along particle 

trajectory is constant (j^ ^ jj). Conservation of magnetic flux (B^dA^ = B^dAj) 

and solid angle area (dA^di’^ - dA^dn^) yields the following relation between initial 

omnidirectional flux (J ) and final omnidirectional flux (J J; from ( J = / j dll ) and 

O * V 47T / 

dA„ Bf \ 

Jf " / jf dUf = / jo dj- dn^ - ^ f j^ dn^ j . 

U ^ ® 4;r / 

J I is obtained. 

During a fnagnetosperic substtrm, many of the plasma injections observed 

are characterized by omnidirectional electron fliOces ner 10^ electrons/ 

cm^-sec. The previous calculations would yield omnidirectional electron 

13 2 

fluxes from high altitude nuclear bursts of approximately Jj. * 10 electrons/cm - 
sec at synchronous orbit. It is clear that these fluxes are significantly greater 
than those encountered during a magnetospheric substorm and would poise a strong 
spacecraft charging ertvironmeht. 
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While the duratlttn of the nuclear burst and resulting injection will be on the 
order of seconds and therefore much shorter than- a typical substorm, studies^® 
indicate the charging process takes only fractions of a second. 


». i ONCLI SIONS 


The atmospheric nuclear burst environment appears to present the potential 
for spacecraft charging effects at synchronous orbits. The calculated fluxes, 
and energies of the Injected electrons are greater than those encountered 
during substorms, and which have been observed to cause spacecraft charging. 

These calculations Contain many approximations, but preliminary results 
indicate that it may be expected that synchronous orbit satellites under certain 
nuclear burst conditions would find themselves subject to a short, but intense, 
period of spacecraft charging. Potentials in the tens of kilovolts are suggested. 
The resulting transient charging and arc discharging Should be a part of the 
design criteria of any spacecraft ekpected to Survive a situation where high 
altitude nuclear bursts are involved. 
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